Transforming a head direction signal into a goal-oriented steering command.

To navigate, we must continuously estimate the direction we are headed in, and we must correct deviations from our goal1. Direction estimation is accomplished by ring attractor networks in the head direction system2,3. However, we do not fully understand how the sense of direction is used to guide action. Drosophila connectome analyses4,5 reveal three cell populations (PFL3R, PFL3L and PFL2) that connect the head direction system to the locomotor system. Here we use imaging, electrophysiology and chemogenetic stimulation during navigation to show how these populations function. Each population receives a shifted copy of the head direction vector, such that their three reference frames are shifted approximately 120° relative to each other. Each cell type then compares its own head direction vector with a common goal vector; specifically, it evaluates the congruence of these vectors via a nonlinear transformation. The output of all three cell populations is then combined to generate locomotor commands. PFL3R cells are recruited when the fly is oriented to the left of its goal, and their activity drives rightward turning; the reverse is true for PFL3L. Meanwhile, PFL2 cells increase steering speed, and are recruited when the fly is oriented far from its goal. PFL2 cells adaptively increase the strength of steering as directional error increases, effectively managing the tradeoff between speed and accuracy. Together, our results show how a map of space in the brain can be combined with an internal goal to generate action commands, via a transformation from world-centric coordinates to body-centric coordinates.

Multimodal cortical neuronal cell type classification.

Since more than a century, neuroscientists have distinguished excitatory (glutamatergic) neurons with long-distance projections from inhibitory (GABAergic) neurons with local projections and established layer-dependent schemes for the ~ 80% excitatory (principal) cells as well as the ~ 20% inhibitory neurons. Whereas, in the early days, mainly morphological criteria were used to define cell types, later supplemented by electrophysiological and neurochemical properties, nowadays. single-cell transcriptomics is the method of choice for cell type classification. Bringing recent insight together, we conclude that despite all established layer- and area-dependent differences, there is a set of reliably identifiable cortical cell types that were named (among others) intratelencephalic (IT), extratelencephalic (ET), and corticothalamic (CT) for the excitatory cells, which altogether comprise ~ 56 transcriptomic cell types (t-types). By the same means, inhibitory neurons were subdivided into parvalbumin (PV), somatostatin (SST), vasoactive intestinal polypeptide (VIP), and "other (i.e. Lamp5/Sncg)" subpopulations, which altogether comprise ~ 60 t-types. The coming years will show which t-types actually translate into "real" cell types that show a common set of multimodal features, including not only transcriptome but also physiology and morphology as well as connectivity and ultimately function. Only with the better knowledge of clear-cut cell types and experimental access to them, we will be able to reveal their specific functions, a task which turned out to be difficult in a part of the brain being so much specialized for cognition as the cerebral cortex.

Whole-brain spatial organization of hippocampal single-neuron projectomes.

Mapping single-neuron projections is essential for understanding brain-wide connectivity and diverse functions of the hippocampus (HIP). Here, we reconstructed 10,100 single-neuron projectomes of mouse HIP and classified 43 projectome subtypes with distinct projection patterns. The number of projection targets and axon-tip distribution depended on the soma location along HIP longitudinal and transverse axes. Many projectome subtypes were enriched in specific HIP subdomains defined by spatial transcriptomic profiles. Furthermore, we delineated comprehensive wiring diagrams for HIP neurons projecting exclusively within the HIP formation (HPF) and for those projecting to both intra- and extra-HPF targets. Bihemispheric projecting neurons generally projected to one pair of homologous targets with ipsilateral preference. These organization principles of single-neuron projectomes provide a structural basis for understanding the function of HIP neurons.

Evolution of neuronal cell classes and types in the vertebrate retina.

The basic plan of the retina is conserved across vertebrates, yet species differ profoundly in their visual needs1. Retinal cell types may have evolved to accommodate these varied needs, but this has not been systematically studied. Here we generated and integrated single-cell transcriptomic atlases of the retina from 17 species: humans, two non-human primates, four rodents, three ungulates, opossum, ferret, tree shrew, a bird, a reptile, a teleost fish and a lamprey. We found high molecular conservation of the six retinal cell classes (photoreceptors, horizontal cells, bipolar cells, amacrine cells, retinal ganglion cells (RGCs) and Müller glia), with transcriptomic variation across species related to evolutionary distance. Major subclasses were also conserved, whereas variation among cell types within classes or subclasses was more pronounced. However, an integrative analysis revealed that numerous cell types are shared across species, based on conserved gene expression programmes that are likely to trace back to an early ancestral vertebrate. The degree of variation among cell types increased from the outer retina (photoreceptors) to the inner retina (RGCs), suggesting that evolution acts preferentially to shape the retinal output. Finally, we identified rodent orthologues of midget RGCs, which comprise more than 80% of RGCs in the human retina, subserve high-acuity vision, and were previously believed to be restricted to primates2. By contrast, the mouse orthologues have large receptive fields and comprise around 2% of mouse RGCs. Projections of both primate and mouse orthologous types are overrepresented in the thalamus, which supplies the primary visual cortex. We suggest that midget RGCs are not primate innovations, but are descendants of evolutionarily ancient types that decreased in size and increased in number as primates evolved, thereby facilitating high visual acuity and increased cortical processing of visual information.

A high-resolution transcriptomic and spatial atlas of cell types in the whole mouse brain.

The mammalian brain consists of millions to billions of cells that are organized into many cell types with specific spatial distribution patterns and structural and functional properties1-3. Here we report a comprehensive and high-resolution transcriptomic and spatial cell-type atlas for the whole adult mouse brain. The cell-type atlas was created by combining a single-cell RNA-sequencing (scRNA-seq) dataset of around 7 million cells profiled (approximately 4.0 million cells passing quality control), and a spatial transcriptomic dataset of approximately 4.3 million cells using multiplexed error-robust fluorescence in situ hybridization (MERFISH). The atlas is hierarchically organized into 4 nested levels of classification: 34 classes, 338 subclasses, 1,201 supertypes and 5,322 clusters. We present an online platform, Allen Brain Cell Atlas, to visualize the mouse whole-brain cell-type atlas along with the single-cell RNA-sequencing and MERFISH datasets. We systematically analysed the neuronal and non-neuronal cell types across the brain and identified a high degree of correspondence between transcriptomic identity and spatial specificity for each cell type. The results reveal unique features of cell-type organization in different brain regions-in particular, a dichotomy between the dorsal and ventral parts of the brain. The dorsal part contains relatively fewer yet highly divergent neuronal types, whereas the ventral part contains more numerous neuronal types that are more closely related to each other. Our study also uncovered extraordinary diversity and heterogeneity in neurotransmitter and neuropeptide expression and co-expression patterns in different cell types. Finally, we found that transcription factors are major determinants of cell-type classification and identified a combinatorial transcription factor code that defines cell types across all parts of the brain. The whole mouse brain transcriptomic and spatial cell-type atlas establishes a benchmark reference atlas and a foundational resource for integrative investigations of cellular and circuit function, development and evolution of the mammalian brain.

Transcriptomic cytoarchitecture reveals principles of human neocortex organization.

Variation in cytoarchitecture is the basis for the histological definition of cortical areas. We used single cell transcriptomics and performed cellular characterization of the human cortex to better understand cortical areal specialization. Single-nucleus RNA-sequencing of 8 areas spanning cortical structural variation showed a highly consistent cellular makeup for 24 cell subclasses. However, proportions of excitatory neuron subclasses varied substantially, likely reflecting differences in connectivity across primary sensorimotor and association cortices. Laminar organization of astrocytes and oligodendrocytes also differed across areas. Primary visual cortex showed characteristic organization with major changes in the excitatory to inhibitory neuron ratio, expansion of layer 4 excitatory neurons, and specialized inhibitory neurons. These results lay the groundwork for a refined cellular and molecular characterization of human cortical cytoarchitecture and areal specialization.

Molecular features driving cellular complexity of human brain evolution.

Human-specific genomic changes contribute to the unique functionalities of the human brain1-5. The cellular heterogeneity of the human brain6,7 and the complex regulation of gene expression highlight the need to characterize human-specific molecular features at cellular resolution. Here we analysed single-nucleus RNA-sequencing and single-nucleus assay for transposase-accessible chromatin with sequencing datasets for human, chimpanzee and rhesus macaque brain tissue from posterior cingulate cortex. We show a human-specific increase of oligodendrocyte progenitor cells and a decrease of mature oligodendrocytes across cortical tissues. Human-specific regulatory changes were accelerated in oligodendrocyte progenitor cells, and we highlight key biological pathways that may be associated with the proportional changes. We also identify human-specific regulatory changes in neuronal subtypes, which reveal human-specific upregulation of FOXP2 in only two of the neuronal subtypes. We additionally identify hundreds of new human accelerated genomic regions associated with human-specific chromatin accessibility changes. Our data also reveal that FOS::JUN and FOX motifs are enriched in the human-specifically accessible chromatin regions of excitatory neuronal subtypes. Together, our results reveal several new mechanisms underlying the evolutionary innovation of human brain at cell-type resolution.

COMBINe enables automated detection and classification of neurons and astrocytes in tissue-cleared mouse brains.

Tissue clearing renders entire organs transparent to accelerate whole-tissue imaging; for example, with light-sheet fluorescence microscopy. Yet, challenges remain in analyzing the large resulting 3D datasets that consist of terabytes of images and information on millions of labeled cells. Previous work has established pipelines for automated analysis of tissue-cleared mouse brains, but the focus there was on single-color channels and/or detection of nuclear localized signals in relatively low-resolution images. Here, we present an automated workflow (COMBINe, Cell detectiOn in Mouse BraIN) to map sparsely labeled neurons and astrocytes in genetically distinct mouse forebrains using mosaic analysis with double markers (MADM). COMBINe blends modules from multiple pipelines with RetinaNet at its core. We quantitatively analyzed the regional and subregional effects of MADM-based deletion of the epidermal growth factor receptor (EGFR) on neuronal and astrocyte populations in the mouse forebrain.

Regulation of early cerebellar development.

The study of cerebellar development has been at the forefront of neuroscience since the pioneering work of Wilhelm His Sr., Santiago Ramón y Cajal and many others since the 19th century. They laid the foundation to identify the circuitry of the cerebellum, already revealing its stereotypic three-layered cortex and discerning several of its neuronal components. Their work was fundamental in the acceptance of the neuron doctrine, which acknowledges the key role of individual neurons in forming the basic units of the nervous system. Increasing evidence shows that the cerebellum performs a variety of homeostatic and higher order neuronal functions beyond the mere control of motor behaviour. Over the last three decades, many studies have revealed the molecular machinery that regulates distinct aspects of cerebellar development, from the establishment of a cerebellar anlage in the posterior brain to the identification of cerebellar neuron diversity at the single cell level. In this review, we focus on summarizing our current knowledge on early cerebellar development with a particular emphasis on the molecular determinants that secure neuron specification and contribute to the diversity of cerebellar neurons.

Study of the potential toxicity of adrenaline to neurons, using the SH-SY5Y human cellular model

Abstract Prolonged overexposure to catecholamines causes toxicity, usually credited to continuous adrenoceptor stimulation, autoxidation, and the formation of reactive pro-oxidant species. Non-differentiated SH-SY5Y cells were used to study the possible contribution of oxidative stress in adrenaline (ADR)-induced neurotoxicity, as a model to predict the toxicity of this catecholamine to peripheral nerves. Cells were exposed to several concentrations of ADR (0.1, 0.25, 0.5 and 1mM) and two cytotoxicity assays [lactate dehydrogenase (LDH) release and 3-(4,5-dimethylthiazol-2yl)-2,5-diphenyl tetrazolium bromide (MTT) reduction] were performed at several time-points (24, 48, and 96h). The cytotoxicity of ADR was concentration- and time-dependent in both assays, since the lowest concentration tested (0.1mM) also caused significant cytotoxicity at 96h. N-acetyl-cysteine (1mM), a precursor of glutathione synthesis, prevented ADR-induced toxicity elicited by 0.5mM and 0.25mM ADR following a 96-h exposure, while the antioxidant Tiron (100µM) was non-protective. In conclusion, ADR led to mitochondrial distress and ultimately cell death in non-differentiated SH-SY5Y cells, possibly because of ADR oxidation products. The involvement of such processes in the catecholamine-induced peripheral neuropathy requires further analysis.

Dual-polarity voltage imaging of the concurrent dynamics of multiple neuron types.

Genetically encoded fluorescent voltage indicators are ideally suited to reveal the millisecond-scale interactions among and between targeted cell populations. However, current indicators lack the requisite sensitivity for in vivo multipopulation imaging. We describe next-generation green and red voltage sensors, Ace-mNeon2 and VARNAM2, and their reverse response-polarity variants pAce and pAceR. Our indicators enable 0.4- to 1-kilohertz voltage recordings from >50 spiking neurons per field of view in awake mice and ~30-minute continuous imaging in flies. Using dual-polarity multiplexed imaging, we uncovered brain state-dependent antagonism between neocortical somatostatin-expressing (SST+) and vasoactive intestinal peptide-expressing (VIP+) interneurons and contributions to hippocampal field potentials from cell ensembles with distinct axonal projections. By combining three mutually compatible indicators, we performed simultaneous triple-population imaging. These approaches will empower investigations of the dynamic interplay between neuronal subclasses at single-spike resolution.

Matrix Inversion and Subset Selection (MISS): A pipeline for mapping of diverse cell types across the murine brain.

The advent of increasingly sophisticated imaging platforms has allowed for the visualization of the murine nervous system at single-cell resolution. However, current experimental approaches have not yet produced whole-brain maps of a comprehensive set of neuronal and nonneuronal types that approaches the cellular diversity of the mammalian cortex. Here, we aim to fill in this gap in knowledge with an open-source computational pipeline, Matrix Inversion and Subset Selection (MISS), that can infer quantitatively validated distributions of diverse collections of neural cell types at 200-µm resolution using a combination of single-cell RNA sequencing (RNAseq) and in situ hybridization datasets. We rigorously demonstrate the accuracy of MISS against literature expectations. Importantly, we show that gene subset selection, a procedure by which we filter out low-information genes prior to performing deconvolution, is a critical preprocessing step that distinguishes MISS from its predecessors and facilitates the production of cell-type maps with significantly higher accuracy. We also show that MISS is generalizable by generating high-quality cell-type maps from a second independently curated single-cell RNAseq dataset. Together, our results illustrate the viability of computational approaches for determining the spatial distributions of a wide variety of cell types from genetic data alone.

Cell-type-specific responses to associative learning in the primary motor cortex.

The primary motor cortex (M1) is known to be a critical site for movement initiation and motor learning. Surprisingly, it has also been shown to possess reward-related activity, presumably to facilitate reward-based learning of new movements. However, whether reward-related signals are represented among different cell types in M1, and whether their response properties change after cue-reward conditioning remains unclear. Here, we performed longitudinal in vivo two-photon Ca2+ imaging to monitor the activity of different neuronal cell types in M1 while mice engaged in a classical conditioning task. Our results demonstrate that most of the major neuronal cell types in M1 showed robust but differential responses to both the conditioned cue stimulus (CS) and reward, and their response properties undergo cell-type-specific modifications after associative learning. PV-INs' responses became more reliable to the CS, while VIP-INs' responses became more reliable to reward. Pyramidal neurons only showed robust responses to novel reward, and they habituated to it after associative learning. Lastly, SOM-INs' responses emerged and became more reliable to both the CS and reward after conditioning. These observations suggest that cue- and reward-related signals are preferentially represented among different neuronal cell types in M1, and the distinct modifications they undergo during associative learning could be essential in triggering different aspects of local circuit reorganization in M1 during reward-based motor skill learning.

Autism genes converge on asynchronous development of shared neuron classes.

Genetic risk for autism spectrum disorder (ASD) is associated with hundreds of genes spanning a wide range of biological functions1-6. The alterations in the human brain resulting from mutations in these genes remain unclear. Furthermore, their phenotypic manifestation varies across individuals7,8. Here we used organoid models of the human cerebral cortex to identify cell-type-specific developmental abnormalities that result from haploinsufficiency in three ASD risk genes-SUV420H1 (also known as KMT5B), ARID1B and CHD8-in multiple cell lines from different donors, using single-cell RNA-sequencing (scRNA-seq) analysis of more than 745,000 cells and proteomic analysis of individual organoids, to identify phenotypic convergence. Each of the three mutations confers asynchronous development of two main cortical neuronal lineages-γ-aminobutyric-acid-releasing (GABAergic) neurons and deep-layer excitatory projection neurons-but acts through largely distinct molecular pathways. Although these phenotypes are consistent across cell lines, their expressivity is influenced by the individual genomic context, in a manner that is dependent on both the risk gene and the developmental defect. Calcium imaging in intact organoids shows that these early-stage developmental changes are followed by abnormal circuit activity. This research uncovers cell-type-specific neurodevelopmental abnormalities that are shared across ASD risk genes and are finely modulated by human genomic context, finding convergence in the neurobiological basis of how different risk genes contribute to ASD pathology.

Virally encoded connectivity transgenic overlay RNA sequencing (VECTORseq) defines projection neurons involved in sensorimotor integration.

Behavior arises from concerted activity throughout the brain. Consequently, a major focus of modern neuroscience is defining the physiology and behavioral roles of projection neurons linking different brain areas. Single-cell RNA sequencing has facilitated these efforts by revealing molecular determinants of cellular physiology and markers that enable genetically targeted perturbations such as optogenetics, but existing methods for sequencing defined projection populations are low throughput, painstaking, and costly. We developed a straightforward, multiplexed approach, virally encoded connectivity transgenic overlay RNA sequencing (VECTORseq). VECTORseq repurposes commercial retrogradely infecting viruses typically used to express functional transgenes (e.g., recombinases and fluorescent proteins) by treating viral transgene mRNA as barcodes within single-cell datasets. VECTORseq is compatible with different viral families, resolves multiple populations with different projection targets in one sequencing run, and identifies cortical and subcortical excitatory and inhibitory projection populations. Our study provides a roadmap for high-throughput identification of neuronal subtypes based on connectivity.

Cell-type specialization is encoded by specific chromatin topologies.

The three-dimensional (3D) structure of chromatin is intrinsically associated with gene regulation and cell function1-3. Methods based on chromatin conformation capture have mapped chromatin structures in neuronal systems such as in vitro differentiated neurons, neurons isolated through fluorescence-activated cell sorting from cortical tissues pooled from different animals and from dissociated whole hippocampi4-6. However, changes in chromatin organization captured by imaging, such as the relocation of Bdnf away from the nuclear periphery after activation7, are invisible with such approaches8. Here we developed immunoGAM, an extension of genome architecture mapping (GAM)2,9, to map 3D chromatin topology genome-wide in specific brain cell types, without tissue disruption, from single animals. GAM is a ligation-free technology that maps genome topology by sequencing the DNA content from thin (about 220 nm) nuclear cryosections. Chromatin interactions are identified from the increased probability of co-segregation of contacting loci across a collection of nuclear slices. ImmunoGAM expands the scope of GAM to enable the selection of specific cell types using low cell numbers (approximately 1,000 cells) within a complex tissue and avoids tissue dissociation2,10. We report cell-type specialized 3D chromatin structures at multiple genomic scales that relate to patterns of gene expression. We discover extensive 'melting' of long genes when they are highly expressed and/or have high chromatin accessibility. The contacts most specific of neuron subtypes contain genes associated with specialized processes, such as addiction and synaptic plasticity, which harbour putative binding sites for neuronal transcription factors within accessible chromatin regions. Moreover, sensory receptor genes are preferentially found in heterochromatic compartments in brain cells, which establish strong contacts across tens of megabases. Our results demonstrate that highly specific chromatin conformations in brain cells are tightly related to gene regulation mechanisms and specialized functions.

Compensatory variability in network parameters enhances memory performance in the Drosophila mushroom body.

Neural circuits use homeostatic compensation to achieve consistent behavior despite variability in underlying intrinsic and network parameters. However, it remains unclear how compensation regulates variability across a population of the same type of neurons within an individual and what computational benefits might result from such compensation. We address these questions in the Drosophila mushroom body, the fly's olfactory memory center. In a computational model, we show that under sparse coding conditions, memory performance is degraded when the mushroom body's principal neurons, Kenyon cells (KCs), vary realistically in key parameters governing their excitability. However, memory performance is rescued while maintaining realistic variability if parameters compensate for each other to equalize KC average activity. Such compensation can be achieved through both activity-dependent and activity-independent mechanisms. Finally, we show that correlations predicted by our model's compensatory mechanisms appear in the Drosophila hemibrain connectome. These findings reveal compensatory variability in the mushroom body and describe its computational benefits for associative memory.

Increased Accuracy to c-Fos-Positive Neuron Counting.

There is not a described method to count the core label of c-Fos-positive neurons, avoiding false-positive and false-negative results. The aim of this manuscript is to provide guidelines for a secure and accurate method to calculate a threshold to select which core of c-Fos-positive neurons marked by immunofluorescence has to be scored. A background percentage was calculated by dividing the intensity value (0 to 255) of the core of c-Fos-positive neurons by its surrounding background from the 8-bit images obtained in a previous study. Using the background percentage from 20% up to 98%, raising 2% once for each score, as threshold to choose which core has to be counted, a script was written for the R program to count the number of the c-Fos-positive neurons and the comparison between control and experimental groups. The differences of the average number of the core counted c-Fos-positive neurons between control and experimental groups, at all thresholds studied, showed a rising value related to an increase of the background percentage threshold as well as a decrease of its p value related to an increase of the threshold of background percentage. For the smallest thresholds (high intensity of label), the differences between groups are suppressed (false negative). However, for the biggest thresholds (nonspecific label), these differences are always the same (false positive). Therefore, to avoid the false-negative and the false-positive values, it was chosen as the threshold of 62% the inflection point of the linear regression, which is equally different from the biggest and smallest values of the differences between groups.

Postsynaptic structure formation of human iPS cell-derived neurons takes longer than presynaptic formation during neural differentiation in vitro.

The generation of mature synaptic structures using neurons differentiated from human-induced pluripotent stem cells (hiPSC-neurons) is expected to be applied to physiological studies of synapses in human cells and to pathological studies of diseases that cause abnormal synaptic function. Although it has been reported that synapses themselves change from an immature to a mature state as neurons mature, there are few reports that clearly show when and how human stem cell-derived neurons change to mature synaptic structures. This study was designed to elucidate the synapse formation process of hiPSC-neurons. We propagated hiPSC-derived neural progenitor cells (hiPSC-NPCs) that expressed localized markers of the ventral hindbrain as neurospheres by dual SMAD inhibition and then differentiated them into hiPSC-neurons in vitro. After 49 days of in vitro differentiation, hiPSC-neurons significantly expressed pre- and postsynaptic markers at both the transcript and protein levels. However, the expression of postsynaptic markers was lower than in normal human or normal rat brain tissues, and immunostaining analysis showed that it was relatively modest and was lower than that of presynaptic markers and that its localization in synaptic structures was insufficient. Neurophysiological analysis using a microelectrode array also revealed that no synaptic activity was generated on hiPSC-neurons at 49 days of differentiation. Analysis of subtype markers by immunostaining revealed that most hiPSC-neurons expressed vesicular glutamate transporter 2 (VGLUT2). The presence or absence of NGF, which is required for the survival of cholinergic neurons, had no effect on their cell fractionation. These results suggest that during the synaptogenesis of hiPSC-neurons, the formation of presynaptic structures is not the only requirement for the formation of postsynaptic structures and that the mRNA expression of postsynaptic markers does not correlate with the formation of their mature structures. Technically, we also confirmed a certain level of robustness and reproducibility of our neuronal differentiation method in a multicenter setting, which will be helpful for future research. Synapse formation with mature postsynaptic structures will remain an interesting issue for stem cell-derived neurons, and the present method can be used to obtain early and stable quality neuronal cultures from hiPSC-NPCs.

DNA methylation atlas of the mouse brain at single-cell resolution.

Mammalian brain cells show remarkable diversity in gene expression, anatomy and function, yet the regulatory DNA landscape underlying this extensive heterogeneity is poorly understood. Here we carry out a comprehensive assessment of the epigenomes of mouse brain cell types by applying single-nucleus DNA methylation sequencing1,2 to profile 103,982 nuclei (including 95,815 neurons and 8,167 non-neuronal cells) from 45 regions of the mouse cortex, hippocampus, striatum, pallidum and olfactory areas. We identified 161 cell clusters with distinct spatial locations and projection targets. We constructed taxonomies of these epigenetic types, annotated with signature genes, regulatory elements and transcription factors. These features indicate the potential regulatory landscape supporting the assignment of putative cell types and reveal repetitive usage of regulators in excitatory and inhibitory cells for determining subtypes. The DNA methylation landscape of excitatory neurons in the cortex and hippocampus varied continuously along spatial gradients. Using this deep dataset, we constructed an artificial neural network model that precisely predicts single neuron cell-type identity and brain area spatial location. Integration of high-resolution DNA methylomes with single-nucleus chromatin accessibility data3 enabled prediction of high-confidence enhancer-gene interactions for all identified cell types, which were subsequently validated by cell-type-specific chromatin conformation capture experiments4. By combining multi-omic datasets (DNA methylation, chromatin contacts, and open chromatin) from single nuclei and annotating the regulatory genome of hundreds of cell types in the mouse brain, our DNA methylation atlas establishes the epigenetic basis for neuronal diversity and spatial organization throughout the mouse cerebrum.